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Abstract 

It is shown that the noise-limited charge sensitivity of a single-electron tran- 
sistor using superconductors (of either SISIS or NISIN type) operating 
near the threshold of quasiparticle tunneling, can be considerably higher than 
that of a similar transistor made of normal metals or semiconductors. The 
reason is that the superconducting energy gap, in contrast to the Coulomb 
blockade, is not smeared by the finite temperature. We discuss also the in- 
crease of the maximum operation temperature due to superconductivity and 
a new peak-like feature on the I — V curve of SISIS structures. 
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Electron transport in the systems of small-capacitance tunnel junctions shows a variety of 
single-electron effects [1]. The simplest and most thoroughly studied circuit revealing these 
effects is the so-called Single Electron Transistor [2] (SET) which consists of two tunnel 
junctions connected in series. At low temperatures (T <C e 2 /Cs, Cs = C\ + C 2 where C\ 
and C 2 are the junction capacitances) the current through this structure depends on the 
background charge Qo of the central electrode (the dependence is periodical with a period 
equal to the electron charge e). Hence, controlling Q Q (for example, by a capacitive gate) it 
is possible to control the current / through the circuit. The possibility to use the SET as a 
highly-sensitive electrometer has been confirmed in numerous experiments. 

The most developed technology of the SET fabrication uses the overlapping narrow alu- 
minum films with a typical junction capacitance about few times 10~ 16 F (see, e.g. Refs. 
[3-5]). Consequently, the operation temperature is typically less than 1 K, and the electrodes 
are in the superconducting state unless the superconductivity is intentionally suppressed by 
the magnetic field. It has been noticed [3,4,6] that the superconductivity of electrodes im- 
proves the performance of the SET (operating near the threshold of quasiparticle tunneling) 
as an electrometer in comparison with the normal-state operation. However, we are not 
aware of quantitative theoretical analysis of this issue, which will be the subject of the 
present paper. 

There are two major characteristics of the SET operation as an electrometer. The first 
one is the amplitude of the output signal modulation for Q variations larger than e. It was 
found experimentally [4] that the use of superconducting electrodes increases the modulation 
amplitude of current / (for fixed bias voltage V) , especially at temperatures comparable to 
e2 /Cs, thus increasing the maximum temperature. The theoretical results of the present 
paper confirm this statement for both NISIN and SISIS structures. 

The other, even more important characteristic of the SET operation is the noise-limited 
sensitivity (ability to detect variations of Q much smaller than e). The best achieved 
sensitivity so far (by the normal state SET) is 7 x I0~ 5 e/ \[Wz at 10Hz [5]. In the present- 
day technology this figure is limited by I/f noise which is most likely caused by random 
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trapping-escape processes in nearby impurities. So, in some sense, the sensitivity does not 
depend much on the parameters of the SET, but rather on the purity of the sample. It is 
unlikely that superconductivity of electrodes can significantly affect these processes. Hence, 
the present-day sensitivities of superconducting and normal SETs with similar parameters 
should not differ much for reasonably low temperatures when both SETs show sufficient 
modulation amplitude. 

However, with the technology improvement one can expect the reduction of the noise due 
to impurities. Then the charge sensitivity of the SET would achieve the limit determined 
by the intrinsic noise [7,8] of the device caused by random electron jumps through tunnel 
junctions (this "white" noise has been recently measured in experiment [9]). Though the 
theory of the "classical" thermal/shot intrinsic noise of the SET is applicable to the gen- 
eral case of one-particle tunneling (normal metals, semiconductors, quasiparticle current in 
superconductors, etc.), most numerical results in Refs. [7] and [8] as well as in a number of 
subsequent papers on this subject (see, e.g. Refs. [10-13]) were obtained only for SETs made 
of normal metals. (Recently some generalization was done [14] to include the possibility of 
two-particle tunneling which can be important in the superconducting case. Let us also 
mention Ref. [6] in which the noise in NISIN SET was briefly considered.) 

In the present paper we apply the theory of Refs. [7] and [8] to the cases of capacitively 
coupled superconducting SI SIS and NISIN SETs (the analysis of a resistively coupled SET 
can be done in a similar way - see Ref. [7]). We show that the noise-limited sensitivity of a 
SET-electrometer can be considerably improved by the use of superconducting electrodes. 

We consider only the quasiparticle tunneling, neglecting the Josephson current, reso- 
nant tunneling of Cooper pairs, Andreev reflection, and cotunneling. This assumption is 
appropriate when the Josephson coupling is negligible and the normal state resistances i?i 
and i? 2 of tunnel junctions are well above the resistance quantum Rq = nh/2e 2 . We use 
the "orthodox" theory [1,2] of the SET and the BCS theory [15] for the calculation of the 
tunneling rates. 

Figure 1 shows the I — V curves at different temperatures for (a) the normal metal 
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NININ case, (b) NISIN case (which is equivalent to SINIS case), and (c)-(d) SI SIS 
case. SETs with C 1 = C 2 and Ri = R 2 = -Re/2 are chosen, and we neglect the gate 
capacitance C g because it can always be formally distributed between C\ and C 2 (see, 
e.g., Ref. [16]). Three curves in each set represent Q = 0, e/4, and e/2, respectively. 
Temperature increase decreases the superconducting energy gap A(T) (which is assumed to 
be equal in all S-electrodes) leading to the noticeable shift to the left of the positions of the 
current jumps in Figs. 1(c) and 1(d). The pure BCS theory would lead to the abrupt jumps 
of the current in SI SIS case. To take into account the unavoidable smoothing of the jumps 
in reality, we assume additionally the inhomogeneous broadening of A(0) with Gaussian 
distribution characterized by the dispersion w . This phenomeno logical parameter is chosen 
as w = 0.05A(0) in Figs. 1(c) and 1(d) (for finite temperatures w(T) = w [A(T)/A(0) - 
(T/A(0))(dA(T)/dT)} was used). 

One can see that in the normal metal case the current / can be considerably modulated 
(ImaJImin £ 2) by Q (V is fixed) only at T < 0.15e 2 /C s , while at T = 0.3e 2 /C s the 
modulation is already negligible, (I max ~ Imin)/Imax — 5%. Notice that the maximum 
relative modulation is achieved at small voltages and does not depend on ratios Ci/C 2 and 
R 1 /R 2 . 

NISIN transistor with A(0) = 0.5e 2 /Cs shows considerable modulation crudely up to 
T « 0.2e 2 /C s , while SISIS transistors with A(0) = 0.5e 2 /C s and A(0) = 2.0e 2 /C E operate 
well almost up to the critical temperature T c (T c /(e 2 /C s ) = 0.28 and 1.14 respectively). The 
case A(0) = 0.5e 2 /Cs corresponds to the typical present-day experimental situation with 
aluminum junctions and C s 0.4 fF (see, e.g. Ref. [4]). Comparison of Figs. 1(c) and 1(d) 
shows that the increase of A(0) provides further improvement of the transistor performance 
at high temperatures. Using Fig. 1(d) one can predict the operation of the niobium-based 
SET with Cs ~ 0.2 fF (current state-of-the-art for aluminum junctions) at temperatures up 
to 7 K. 

Superconductivity improves the SET performance at relatively high temperatures be- 
cause, in contrast to the Coulomb blockade, the superconducting energy gap is not smeared 
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by the finite temperature. In the normal metal case the I — V curve has a cusp at the 
Coulomb blockade threshold V t = mm. i>n {V^ n \ V i>n > 0}, where 



and this cusp is rounded within the voltage interval proportional to the temperature. In 
SISIS case the jump of the I — V curve at V t which is shifted due to the energy gap, 
V t = mm iin {V iin + 2A(T)C E /eCi | V t > 4A(T)}, remains sharp even at T ~ A(T), 
and the subthreshold current increase is only proportional to exp(— T/A(T)). This ex- 
plains why SISIS transistor shows considerable dependence on Q for the temperatures 
almost up to T c even if T > e 2 /C s . In NISIN case the I-V curve in the vicinity of 
V t = min ijri {Vi jn + A(T)C s /eCj | V t > 2A(T)} is rounded by the finite temperature, that 
makes NISIN transistor worse than SISIS transistor, however, it is still better than usual 
NININ transistor. 

Let us briefly discuss the origin of small peaks of the current at moderate temperatures 
visible in Figs. 1(c) and 1(d) (SISIS case) at voltages close to the middle of the subthreshold 
region. The position of a peak satisfy Eq. (1) and corresponds to zero energy gain W for a 
particular tunneling process (hence, it coincides with the position of one of the I — V cusps in 
the corresponding NININ SET). In this case the singularities in the density of states of two 
electrodes match, leading to increase of tunneling of thermally excited quasiparticles. Hence, 
the origin of peaks is similar to that of well-known peaks [15] at V — (A 1 (T) — A 2 (T))/e 
in the single junction made of superconductors with different gaps A X (T) and A 2 (T). In 
our case energy gaps are the same but the Coulomb blockade provides the relative shift of 
the singularities in the density of states. The analysis of the master equation [1,2] shows 
that the singularity-matching peak can be significantly high only within the voltage range 
2A(T) < V < 2A(T) + e/C s . Hence, not more than two closely located peaks from the set 
(1) can be well pronounced on the I — V curve. In the case A(T) > e 2 /Cs the peak position 
is close to the center of the subthreshold region, and, hence, close to the peak due to the 
Josephson-plus-quasiparticle cycle [17]. 




(1) 
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Peculiarities of another type (jumps of the current) also exist at V < V t at finite tem- 
peratures. Their positions satisfy equation V = 2A(T)Cs/eCj + V ijn and correspond to 
the energy gain W = 2A(T) for a particular tunneling process. The jump height decreases 
with the decrease of voltage and vanishes at V < 2A(T)/e. These jumps of current are not 
well- noticeable in Figs. 1(c) and 1(d). 

Now let us consider the noise-limited sensitivity of the SET. The minimum detectable 
charge for the given bandwidth Af is 5Q = (S/A/) 1 / 2 / (dl/dQo) where the spectral density 
Si of the current noise is taken in the low frequency limit. The ultimate low-temperature 
(T < e 2 /C s ) sensitivity in the NININ case is [7,8] min5Q ^ 2.7Cx{R min TAf) 1 / 2 , R min = 
mm{R 1: R 2 }. This result can be somewhat improved in the NISIN SET (with the same 
resistances) operating near the threshold V t of quasiparticle tunneling. At low temperatures, 
T <C min{e 2 /Cs, A(T)}, and for V close to nondegenerate V t , we can use approximation 
Si ~ 2eJ, / ~ h^V -Vt^C^dC^), where I^v) = (1/ eR i )[TA(T)/2] 1 / 2 / °° dyj v /y/[l + 
exp(y + (A — ev)/^)}" 1 is the "seed" I — V curve of i-ih junction. Then the ultimate 
sensitivity is given by equation 

min<5Q = C s (2eA/) 1 / 2 mm{y/l^v)/(dl /dv)}, (2) 

and finally we get the result 

min5g ~ 2.6C s (i? mm TA/) 1 / 2 [T/A(T)] 1 / 4 (3) 

which is better than NININ sensitivity when T < A(T). The main reason for the improve- 
ment is the increase [3,4,6] of the transfer coefficient dI/dQ ~ (e/Ci)(dI /dV), because the 
differential resistance Rd of the "seed" I — V curve near the onset of quasiparticle tunneling 
is less than Ri. Notice that the "orthodox" theory used here is valid only if R d > Rq because 
the cotunneling processes [18,6] impose the lower bound for (dl / dV)" 1 on the order of Rq 
[19]. For relatively high temperatures the ratio of minimum SQ in NISIN and NININ 
cases is larger than [A(T)/T] 1//4 (e.g., compare the dashed lines in Fig. 2) because NININ 
sensitivity starts to deviate up from the low-temperature approximation at smaller T than 
NISIN sensitivity. 
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The improvement of the ultimate sensitivity is more significant in SI SIS SET. For pure 
BCS model the "orthodox" theory gives infinite derivative dI/dQ at V — V t even for finite 
temperature leading to 5Qo — > 0. Hence, the "orthodox" ultimate sensitivity depends on the 
imperfection of the current jump which is described in our model by the energy gap spread 
w ( Wo <min{A(T),e 2 /C s }). 

Figure 2 shows 5Qo together with current / and ratio Si/2el, as functions of the voltage 
for the symmetric SISIS SET with parameters A(0) = 0.5e 2 /C s , w = 0.05A(0), T = 
0.1e 2 /Cs, and Q = 0.25e (numerical calculations are done using the method described in 
Refs. [7] and [8]). Dashed lines show 5Qo for similar NININ and NISIN SETs. One can 
see that the sensitivity of SISIS SET is much better than for NININ and NISIN cases 
within a relatively narrow voltage range which corresponds to the jump of current. 

In contrast to NININ and NISIN cases, the approximation Si ~ 2el is not accurate 
in the vicinity of V t for SISIS SET even at low temperatures (see Fig. 2) because the 
relatively large tunneling rate in the junction determining V t , is comparable to the tunneling 
rate in the other junction. This approximation is valid only if T <C A(T) <C e 2 /Cs, and 
would lead to inaccuracy typically about 10% for the analytical calculation of min^Qo if 
T <C A(T) ~ e 2 /C s . Nevertheless, it can be used as a crude estimate. Using Eq. (2) and 
smoothed by wo low-temperature (T <C A(T)) "seed" I — V curve for SIS junction [15] we 
get 



Notice that the numerical factor depends on the particular model describing the shape of 
the current jump. Comparing Eq. (4) with the result for NININ SET, we see that the 
temperature T is replaced in SISIS case by Wq/A(T). Hence, the ultimate sensitivity is 
better in SISIS SET (resistances are the same) with sufficiently narrow width of the current 



In the case of very sharp "seed" I — V curve, w < A(T)Rq/ R iy the slope of the jump of 
the SET I — V curve is determined by cotunneling [18] and it cannot be sharper than crudely 




(4) 



jump, w < (TA(T))V2. 
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Rq 1 [19]. Then mmSQ is on the order of C S (A/ A(T)R 2 Q / R) 1 / 2 (we assume A(T) > e 2 /C s , 
R 1 = R 2 ), and the ultimate sensitivity is better than for NININ SET if T > A{T){R Q /R) 2 . 
The sensitivity of such an ideal SISIS SET is even better than the "quantum" (T = 0) 
sensitivity of a symmetric = i? 2 ) NININ SET operating at Vf ~ e/Cs (in that case [7] 
min^Q ~ (hC^ A/) 1//2 ), if R/Rq ^ A(T)C s /e 2 . However, notice that the quantum-noise- 
limited min^Q of a NININ SET can be made arbitrary small using small V t (and large 
ratio R/Rq) [7] or large ratio R1/R2 [8]; hence, in this sense the use of superconducting 
electrodes cannot improve further the ultimate sensitivity. 

The author thanks D. V. Averin and K. K. Likharev for valuable discussions and critical 
reading of the manuscript. The work was supported in part by ONR grant N00014-93-1-0880 
and AFOSR grant 91-0445. 
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FIGURES 

FIG. 1. I - V curves for (a) NININ, (b) NISIN (or SINIS), and (c)-(d) SISIS SETs 
for three values of Qq (0, e/4, and e/4) and several temperatures T. The curves for different T 
are offset vertically for clarity. Notice that the modulation by Qo survives up to higher T in the 
superconducting transistors. 

FIG. 2. The minimum detectable charge SQo, the current /, and the ratio Sj/2el as functions 
of the bias voltage V for SISIS SET. Dashed lines show 5Q for NININ and NISIN SETs. The 
best sensitivity is achieved in SISIS case. 
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Figs, la, b 
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